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Abstract 

We present a collection of matrix valued shape invariant potentials which 
give rise to new exactly solvable problems of SUSY quantum mechanics. It 
includes all irreducible matrix superpotentials of the generic form W = kQ + 
j:R + P where k is a variable parameter, Q is the unit matrbc multiplied by a 
real valued function of independent variable x, and P, R are hermitian matri- 
ces depending on x. In particular we recover the Pron'ko-Stroganov "matrbc 
Coulomb potential" and all known scalar shape invariant potentials of SUSY 
quantum mechanics. 

In addition, five new shape invariant potentials are presented. Three of 
them admit a dual shape invariance, i. e., the related hamiltonians can be fac- 
torized using two non-equivalent superpotentials. We find discrete spectrum 
and eigenvectors for the corresponding Schrodinger equations and prove that 
these eigenvectors are normalizable. 
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1 Introduction 



Invented by E. Witten [JJ as a toy model supersymmetric quantum mechanics (SSQM) 
became a fundamental field including many interesting external and internal prob- 
lems. In particular the SSQM presents powerful tools for explicit solution of quan- 
tum mechanical problems using the shape invariance approach [2]. Unfortunately, 
the number of problems satisfying the shape invariance condition is rather restricted. 
However, such problems include practically all cases when the related Schrodinger 
equation is exactly solvable and has an explicitly presentable potential. Well known 
exceptions are exactly solvable Schrodinger equations with Natanzon potentials [5] 
which are formulated in terms of implicit functions. The list of shape invariant po- 
tentials depending on one variable can be found in |Q. 

An interesting example of QM problem which admits a shape invariant supersym- 
metric formulation was discovered by Pron'ko and Stroganov [3] who studied a motion 
of a neutral non-relativistic fermion which interacts anomalously with the magnetic 
field generated by a thin current carrying wire. 

The supersymmetric approach to the Pron'ko-Stroganov (PS) problem was first 
applied in paper [6] with using the momentum representation. In paper [7] this prob- 
lem was solved using its shape invariance in the coordinate representation. Recently 
a relativistic generalization of the PS problem was proposed [S] which can also be 
integrated using its supersymmetry with shape invariance. 

The specificity of the P S problem is that it is formulated using a matrix super- 
potential while in the Standard SSQM the superpotential is simply a scalar function. 
Matrix superpotentials themselves were discussed in many papers, see, e. g., [9], [TU] . 
[TT] . [T3] but this discussion was actually reduced to analysis of particular examples. 
In papers [U] such superpotentials were used for analysis of motion of a spin \ par- 
ticle in superposed magnetic and scalar fields. In paper [llj a certain class of such 
superpotentials was described which however was ad hoc restricted to 2 x 2 matrices 
which depend linearly on the variable parameter. Thus, in contrast to the case of 
scalar superpotentials, the class of known matrix potentials includes only few exam- 
ples which are important but rather particular, while the remaining part of this class 
is still "terra incognita". It seems to be interesting to extend our knowledge of these 
potentials since this way it is possible to find new systems of Schrodinger equations 
which are exactly integrable. 

In the present paper a certain class of matrix valued superpotentials is described 
which includes the superpotential for the PS problem as a particular case. Moreover, 
we do not make a priori suppositions about the dimension of the involved matri- 
ces but restrict ourselves to linear and inverse dependence of the superpotentials on 
variable parameters. Trough our approach however the problem of classification of 
indecomposable matrix potentials which are shape invariant appears to be completely 
solvable. We present a classification of such potentials and discuss the corresponding 
exactly solvable problems for coupled systems of Schrodinger equations. In partic- 
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ular the discrete energy spectra and exact solutions for these models are found and 
normalizability of the ground and exited states is proven. Solutions corresponding to 
continuous spectra are not considered here. 

Three out of flve found hamiltonians admit alternative factorizations with using 
different superpotentials. The corresponding potentials are shape invariant w. r. t. 
shifts of two different parameters. Such dual shape invariance results in existence of 
two alternative spectra branches. Moreover, for some values of free parameters both 
these branches can be realized. 



2 Superpotential for PS problem 

The PS problem was discussed in numerous papers, see, e. g., [5]- [7J. Thus we will 
not present its physical motivations and calculation details but start with the corre- 
sponding equation for radial functions [7J 

H K ip = E k ^j (1) 

where H K is a Hamiltonian with a matrbc potential, E K and ip are its eigenvalue 
and eigenfunction correspondingly, moreover, ip is a two-component spinor. Up to 
normalization of the radial variable x the Hamiltonian H K can be represented as 

d 2 1 1 

= -Qtf + < K - a ^~i + ^~ ( 2 ) 

where a\ and 03 are Pauli matrices and k is a natural number. In addition, solutions 
of equation ([T]) must be normalizable and vanish at x = 0. 
Hamiltonian H K can be factorized as 

= a+a~ + c K (3) 

where 

a *=iL +w - a «=-iL +w - ck= -(2tti7 

and W is a matrix superpotential 

1 1 2/c + l 

W K = —a 3 - —— -o-i - — • 4 

2x 2k + 1 2x 

Another nice property of Hamiltonian H K is that its superpartner is equal to 
H K+ x, namely 

d 2 11^ 

H K = a K a K + C K = + ( K + 1 )( K+ 1 ~ °3)^ + <Tl- = H K + l 

Thus equation ([TJ admits supersymmetry with shape invariance and can be solved 
using the Standard technique of SSQM [3]. 
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3 Generic matrix shape invariant potentials 

Following a natural desire to find other matrix potentials which are form invariant 
we consider equation ([T]) with 

H k = -^ + V k (x), (5) 

where V k (x) is an n x n dimensional matrix potential depending on variable x and 
parameter k. 

Suppose that the Hamiltonian accepts factorization fl3]) where W k (x) is a super- 
potential. Our goal is to find such superpotentials which generate form invariant 
potentials V k (x). 

Assume W k (x) is Hermitian. Then the corresponding potential V k (x) and its 
superpartner V k + (x), i. e., 

V k (x) = -^ + W> and V+(x) = ^ + WZ (6) 

are Hermitian too. 

Suppose also that the Hamiltonian be shape invariant, i. e., 

H k = H k+a + C k , (7) 

thus V k + = V k+a + C k or 

Wl + W' k = W 2 k+a - W' k+a + C k (8) 

were C k and a are constants. 

Let us state the problem of classification of shape invariant superpotentials, i. e., 
n x n matrices whose elements are functions of x, k satisfying conditions ([S]). In the 
following section we present such classification for a special class of superpotentials 
whose dependence on k is defined by terms proportional to k and r only. 

4 Irreducible matrix superpotentials 

To generalize (j3J) we consider superpotentials of the following special form 

W k = kQ + ^R + P, (9) 

where P, R and Q are n x n Hermitian matrices depending on x. Moreover, we 
suppose that Q = Q(x) is proportional to the unit matrix. This supposition can be 
motivated by two reasons: 
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• our goal is to generalize superpotential fll]) in which the term linear in k is 
proportional to the unit matrix; 

• restricting ourselves to such Q it is possible to make a complete classification of 
the corresponding superpotentials satisfying shape invariance condition (jHJ). 

We do not make any a priori supposition about possible values of the continuous 
independent variable x. However we suppose that relations ([8]) are valid also for 
k — > k' where k 1 = k + a, k + 2a, k + na and n is a natural number which is 
sufficiently large to make the following speculations. 

It is reasonable to restrict ourselves to the case when the matrices P and R 
cannot be simultaneously transformed to a block diagonal form since if such (unitary) 
transformation is admissible, the related superpotentials are completely reducible. 
Thus we suppose that the pair of matrices < P, R > is irreducible. Let us show that 
in this case it is sufficient to consider 1x1 and 2x2 matrices only. 

Considering the special case a = we conclude that the corresponding W k should 
be linear in x provided relation ([8]) is satisfied: 



where C k is a constant multiplied by the unit matrix and M k is a constant hermitian 
matrix which can be diagonalized. In this way we obtain a direct sum of shifted one 
dimensional oscillators whose irreducible components can be represented in form (Ti~0]) 
where C k and M k are constant s. 

Let a ^ 0. Substituting ([9]) into ([8]) and multiplying the obtained expression by 
k 2 (k + a) 2 we obtain 



AB 2 (Q' - aQ 2 ) + 2B 2 (P' - aQP) + aB{R, P} + ABR' + aAR 2 = B 2 C k (11) 



where {R, P} = RP + PR is anticommutator of matrices P and Q, A = 2k + a, 
B = kik + a) and the prime denotes derivative w. r. t. x. 

Ali terms in the l.h.s. of equations ffTT]) are polynomials in discrete variable k. In 
order for this equation be consistent, its r. h. s. (which includes an arbitrary element 
Ck) should also be a polynomial of the same order whose general form is 



where the Greek letters denote arbitrary parameters. Substituting (Fl2l into (ITT]) and 
equating coefficients for linearly independent terms we obtain the following system: 



W k = -C k x + M k 



(10) 



B 2 C k = volAB 2 - 2fiB 2 - a\B + pAB + au 2 A 



(12) 



Q 2 a - Q' + va = 0, 
P' - aQP + /i = 0, 
{R,P} + \ = 



(13) 
(14) 
(15) 
(16) 



R' = p, R 2 = co 2 . 
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It follows from ffTB]) that p = and R is a constant matrix whose square is 
proportional to the unit one. 

If R is proportional to the unit matrix J or is the zero matrix (in the last case 
u = 0) the corresponding superpotential (jHJ) is reducible. Let u 7^ and i? ^ ±ul 
then the general form of P satisfying (I15p is 

P = A i? + J p (17) 

where P is a matrix which anticommutes with R. 

A straightforward analysis of equation (j!4j) shows that it is easily integrable, but 
to obtain non-trivial Q it is necessary to set /i = A = 0. Indeed, without loss of 
generality hermitian matrix R whose square is proportional to the unit matrix can 
be chosen in the diagonal form: 

r, f Imxm mX s l | /1 o\ 

R = lo \ T I, m + s = n (18) 

where u ^ is a constant, I and 0... are the unit and zero matrices whose dimension 
is indicated in subindices, and without loss of generality we suppose that s > m. 
The corresponding matrbc P satisfying fliZo'j) has the following generic form: 

\ M t J { ' 

\ lvl sxm u sxs / 

where M mxs is an arbitrary matrbc of dimension m x s. Substituting (fT7j) - (|T9ll into 
fTl4]) we obtain the following equations: 

P' = 2aQP, (20) 

X R + fiI nxn ) Q = 0. (21) 



Analyzing equation f l2T|) we conclude that for A 2 + /i 2 7^ the matrbc in brackets 
is invertible and so we have to set Q = 0. If Q is nontrivial we have to set A = fi = 0. 
As a result the system ( fl4]) - (fr6|) is reduced to the following form 

Q 2 a - Q' + vet = 0, (22) 

P = P exp (a IQdx) , (23) 



{P,P} = 0, R z = u\ (24) 

_ aa; 2 (2fc + a) . , , . 

Ct= F (fc + a p +TO(2<: + a) (25) 

where both P and P are constant matrices. 
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Thus the problem of classification of matrix valued shape invariant potentials ([H]) 
is reduced to solving the first order differential equation (|22j) for function Q and the 
algebraic problem (124]) for hermitian matrices R and P. 

Let us show that hermitian n x n matrices P and R which satisfy conditions (l24j) 
can be simultaneously transformed to a block diagonal form. Moreover, irreducible 
matrices satisfying ( 12"4"|) are nothing but the 2x2 Pauli matrices multiplied by con- 
stants, and "1x1 matrices" (scalars) satisfying RP = 0. Starting with (TTB|) and (JT^J) 
and applying a unitary transformation 



R^R' = URU\ P ->■ P' = PPP f , 
P = 



r U>mxm m xs 
0«vm Wovo 



p , = , x W M^ xs R' — R (26) 



where u mxm and u sxs are unitary submatrices, we obtain 

0mxm -A 

/if't D 

iw sxrra u sxs 

with 

M mxs — u mxm M mXs u sxs . (27) 

Transformation (J27j) can be used to simplify submatrix M mxs . In particular this 
submatrix can be reduced to the following form 



^Cxa = ^mxm mx (s-m)J (28) 

where M mxm is a diagonal matrbc: 

M mxm = diag(//i,/^ 2 , • • • ,/i m ) (29) 

where /z 2 , ... are real parameters. Without loss of generality we suppose that there 
are r nonzero parameters yUi,yU2, ■••AV with < r < m being the rank of matrbc M. 

Notice that transformation f[2"Tl) - f[2"9l for rectangular matrices M is called singular 
value decomposition. Such transformations are widely used in linear algebra, see, e. g., 

PU- 

But the set of matrices {R, P A } with R and P A given in (fT8"D and ((2SD, (12gj). (1291 
is completely reducible since by an accordant permutation of rows and columns they 
can be transformed to direct sums of 2 x 2 matrices {P 2 x2, P2x2} where 

P2X2 = ( q ^ u J=U<T 3 , P2X2 = ^ Q^=/iCTi, /i = /ii,/i 2 , (30) 

and of 1 x 1 matrices 

Pixi = ±cj, Axi = /i, nu = (31) 
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were u and /i are arbitrary real numbers. The transformation of matrices ( TTSj) and 
( l2"Ej) . ( I2"gj) . ( I2"9"j) to the direct sum of matrices ( 13T)j) and (13"Tj) can be given explicitly as 

i? ->• PPP f , P ->• C/PC/T 

where C/ is a unitary matrix whose nonzero entries are: 

U a a = Ub m+b-l = U m+ b 6+1 = 1; 

a = 1, m + s, m + s + 1, m + s + 2, • • • , n, b = 2, 3, • • • , s + 1. 

Thus up to unitary equivalence we have only two versions of irreducible matrices 
P and P which are given by equations (13"U1) and (1311) . 

The remaining equation (1221) is easily integrable, thus we can find all inequivalent 
irreducible superpotentials (Q in explicit form. 



5 Superpotentials and shape invariant potentials 

There are six different types of solutions of equation (T2"2"l) . namely 

Q = 0, v = 0, (32) 

and 



Q = 


(XX 






Q = 


X X 2 

a a 2 


<0, 




Q = 


A , ^ 

tanh Xx, v 

a 


A 2 
a 2 


<0, 


Q = 


coth Xx, v 

a 


X 2 

a 2 


<0, 


Q = 


X 

— tan Xx, v = 
a 


X 2 

— >0 





that are defined up to translations x — > x + c, c is an integration constant, and a is 
supposed to be nonzero. 

The corresponding matrices P are easily calculated using equations ( 12"3"j) and (150)1 
or (13T]1. 

Let us note that using solutions (l32|) or solutions ( )33|) for scalar P and P given 
by relations Q31]L we simply recover the known list of shape invariant potentials 
which is presented, e. g., in [J], see the table on pages 291-292 (this list includes also 
the harmonic oscillator ( TIU]) ). We will not present this list here but note that our 
approach gives a simple and straightforward way to find it. 

Consider the case when P and P are 2x2 matrices ( |30|) . Now solutions (132]) for 
Q and solutions with trivial matrices P are not available since they lead to reducible 
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superpotentials. However, solutions f[3"3"j) are consistent. Substituting f l2~3"j) . ( )3"0"]) . 
into ([H]) we obtain the following list of matrix superpotentials 

W^ = ((2 f i + l)a 3 -2K-l)j- + ^-^a 1 , M > ~, (34) 

W K)IX = A (-K + fiexp(-Xx)a 1 - ^cr 3 ) , (35) 

W Kj h = A (k tan Aa; + /i sec Axcr 3 + — <7iJ , (36) 

W/c,/* = A ^— KcothAx + /zcschAxcr 3 — — cr^ , /i < 0, a; > 0, (37) 

Wk,h — A (—k tanh Ax + /i sech Axo"i er 3 j , (38) 

where we introduce the normalized parameter k = -. These superpotentials are de- 
fined up to translations x — )■ x + c, /t — > k + 7, and up to unitary transformations 
W K ^ — >■ U a W K ^Ul where £/1 = o"i, C/2 = 775(1 =t lcr 2) and f/ 3 = <t 3 . In particular these 
transformations change signs of parameters \i and u in fl35]) -f l38]) and of /i + \ in fl34|) , 
thus without loss of generality we can set 

u > 0, /i > (39) 

in all superpotentials fl35|) - fl38|) . Zero values of these parameters are excluded if 
superpotentials (j3i|) - (j3T]l are irreducible. 

Conditions (I3"9"j) can be imposed also for superpotential (J3"7j) . To unify some fol- 
lowing calculations we prefer to fix the signs of /i and n in the way indicated in 

fl3U). 

Notice that the transformations k — > k' = k + a correspond to the following 
transformations for k: 

K ^ K ' = K + l. (40) 

If /i = and u — 1 then operator (1341 coincides with the well known superpoten- 
tial for PS problem (TJJ, but for /i 7^ superpotential ( 1541) is not equivalent to fll]). The 
other found superpotentials are new also and make it possible to formulate consis- 
tent, exactly solvable problems for Schrodinger equation with matrbc potential. The 
corresponding potentials V K can be found starting with fl34|) - fl37|) and using defmition 
([H]). Let us rewrite equation ([6]) as follows: 

K, ~ K, = V K = V K + c K (41) 

where c K is a constant and V K does not include constant terms proportional to the 
unit matrbc. As a result we obtain 

V K = (n(n + 1) + k 2 - k{2h + l)a 3 ) — - -a u (42) 
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V K = \ 2 exp(-2Ax) - (2k - exp(-Ax)tr 1 + 2wct 3 ) , (43) 
V« = A 2 [(k(k — 1) + fJ 2 ) sec 2 \x + 2co tan \xa\ 
+/i(2k — 1) sec Ax tan Axcr 3 ) , 
Vk = A 2 — 1) + /i 2 ) csch 2 (Ax) + 2wcothAa;cri 

(45) 

— 2k) coth Ax csch Axo"3) , 
Vk = A 2 ((;U 2 — — 1)) sech 2 Ax + 2wtanh Axo"3 

(46) 

— /x(2k — 1) sech Axtanh Axai) . 

Potentials (H2]l . fHB jl . ( H3)l (g5]) and (06]) are generated by superpotentials (jSD , fl35|) . 
( I36|) . fl37|) and fl38|) respectively. The corresponding constants c K in f )4T|) are 



u; 2 



(2k + 1) 2 
for potential (142]) . 



(47) 



A 2 (. 2 + ^) (4* 



for potentials fl43|) . (]45|). f!46|) and 



C = A 2 - A (49) 
for potential fT44|) . 

Ali the above potentials are shape invariant and give rise to exactly solvable prob- 
lems for systems of two coupled Schrodinger equations, i. e., for systems of Schrodinger- 
Pauli type. 



6 Dual shape invariance 

To find potentials (|42|) -( 145|) we ask for their shape invariance w. r. t. shifts of parameter 
k. The shape invariance condition together with the supposition concerning the 
generic form of the corresponding superpotential make it possible to define these 
potentials up to arbitrary parameters A, u, k and /z. 

Starting with superpotentials fl34l) - fl37|) we can find the related potentials f )42|) - 
f )45|) in a unique fashion. But let us consider the inverse problem: to find possible 
superpotentials corresponding to given potentials wich in our case are given by equa- 
tions fT42l) - fr45l) . The problems of this kind are very interesting since their solutions 
can be used to generate families of isospectral hamiltonians. It happens that in the 
case of matrbc superpotentials everything is much more interesting since there exist 
additional superpotentials compatible with the shape invariance condition. 
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To find the mentioned additional superpotentials we use the following observation: 
potentials fj42H . fj44jl and fj45H are invariant with respect to the simultaneous change 

11 

-►«--, (50) 

In addition, there exist another transformations of /i and k but they lead to the same 
results. 

Thus in addition to the shape invariance w. r. t. shifts of k potentials (|42p . (|44p 
and fT4"5]) should be shape invariant w. r. t. shifts of parameter fj, also. In other words, 
superpotentials in Section 5, should be considered together with superpotentials which 
can be obtained from (IMjl . ( 1361) and (|37|) using the change (150]) . 

Thus, we also can represent potentials (I34p . (1361) and (1371) in the following form 

K K - W'^ = V, + c, (51) 

where = V K , and 

~ Kcr 3 — /i — 1 uj uj 1 

W ^ = x + WTT)^ C » = W^? (52) 

for Vfc given by equation ( l4"2"|) . 

VV MjK = ^ ^(2/i + 1) tan Xx + (2k - 1) sec Xxa 3 + 2 ^ (53) 
for potential (T4"4"p . and 

W^ K = ^ f-(2/i + 1) coth Ax + (2k - 1) csch Ax<x 3 - ^^J a ^j ( 54 ) 
for potential (T4"5"]) . The related constant constant is: 

^ = ^(±^+1)' + ^^) (55) 

where the sign " +" and " — " corresponds to the cases flo^j) and (IMI) respectively. 

We stress that superpartners of potentials (1511) constructed using superpotentials 
W^, i.e., 

V; = W*, + W'^ (56) 
satisfy the shape invariance condition since 

V+ = K+i + C, 
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with = c^+i - c M . 

Thus potentials f )34|) . fl36l) and fl37j) admit a dual supersymmetry, i. e., they are 
shape invariant w. r. t. shifts of two parameters, namely, k and /i. More exactly, 
superpartners for potentials f H2|) . fj44|) and f H5j) can be obtained either by shifts of 
k or by shifts of fi while simultaneous shifts are forbidden. We call this phenomena 
dual shape invariance. 

Notice that the remaining potentials (|43|) and ( H6|) do not posses the dual shape in- 
variance in the sense formulated above. In potential (H3|) parameter /i is not essential. 
It is supposed to be non-vanishing (since for ji = the corresponding superpotential 
is reducible) and can be normalized to the unity by shifting independent variable x. 

The hamiltonian with potential f )46p is not invariant w. r. t. change f )50p . However 
if we suppose that parameter /i be purely imaginary, i. e., set /i = i/i with ji real, the 
corresponding potential admits discrete symmetry fl50|) for parameters k and ji and 
thus possesses dual supersymmetry with shape invariance. In this way we obtain a 
consistent model of " PT-symmetric quantum mechanics [15]" with the dual shape 
invariance. Discussion of this model lies out the scope of present paper. We only note 
that for u = the corresponding potential is decoupled to a direct sum of potentials 
discussed in [16]. 



7 Exactly solvable problems 

for systems of Schrodinger equations 

Consider the Schrodinger equations 

= (-^ + V*) V = E K tfj (57) 

where H K = o^^a'^ + c,,, and V K are matrix potentials represented in (H5|) - (H5|) . Since 
all these potentials are shape invariant, equations (I57p can be integrated using the 
Standard technique of SSQM. An algorithm for construction of exact solutions of 
supersymmetric a shape invariant Schrodinger equations includes the following steps 
(see, e.g., [1]): 

• To find the ground state solutions ipo( 

) which are proportional to square 
integrable solutions of the first order equation 

d_ 

dx 



a K,^o(K, fi,x)= + W Ktfl ) 1p (K, /i, X) = 0. (5J 



In view of f p£Tj) function ip (n,fi,x) solves equation (JSTjl with 

E K = E Kfi = -c K . (59) 



11 



• To find a solution ipi( ) for the first excited state which is defined by the 
following relation: 

V>i(k, /i, x) = o^Vo(« + 1, A 1 ; a;) = (~ J~ + ^«.^ ^oO + 1> /*> f)- (60) 

Since and satisfy the interwining relations 

H* a K,n = a K4l H K+1 
function (EDJ) solves equation ( 157)) with = = — c K+ i. 

• Solutions for the second excited state can be found as ip2(K,/JL,x) = a+ + 
l,fi,x), etc. Finally, solutions which correspond to n th exited state for any 
admissible natural number n > can be represented as 

ip n {K, fi, x) = a^a+ +1>/i • • • a+ +n _ x ^ {k + n, fi, x). (61) 

The corresponding eigenvalue E Rjn is equal to — c K+n . 

• For systems admitting the dual shape invariance it is necessary to repeat the 
steps enumerated above using alternative (or additional) superpotentials. 

Ali potentials presented in the previous section generate integrable models with 
Hamiltonian (157)) . However, it is desirable to analyze their consistency. In particular, 
it is necessary to verify that there exist square integrable solutions of equation (158)) 
for the ground state. 

In the following sections we prove that such solutions exist for all superpotentials 
given by equations ( 153)) -( 157)) and ( 152)) - ( 153)) . We will see that to obtain normalizable 
ground state solutions it is necessary to impose certain conditions on parameters of 
these superpotentials. 

To finish this section we present energy spectra for models (157)) with potentials 
(S2])-g5]): 



for potential (132)) . 



E = -X'{N* + £) (63) 



N 2 J 

for potentials (133")) , (13B)) . (135)) . and 

E = x ' ( N2 ~ w) (64) 
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for potentials (jUJ). 

In equations f[6"2"j) - f[6"2"j) we omit subindices labeling the energy levels. The spectral 
parameter N can take the following values 

N = n + n, (65) 

and (or) 

N = n + /i + - (66) 

where n = 0,1,2,... are natural numbers which can take any values for potentials 
P2]) - p3]l . For potentials f}4"3"]h fl46|) and (@5D with a fixed k < the admissible values 
of n are bound by the condition (k + n) 2 > see section 9. 

For potential (j4"3"j) the spectral parameter is defined by equation fl65|) . For poten- 
tials (T4"2"j) . (|44|) . (|46p the form of ./V depends on relations between parameters k and 
/i, see section 9. 



8 Some special values of parameters and isospec- 
trality 

Let us show that for some values of parameters /i and k potentials Q42p - ([4"6"]) are 
isospectral with direct sums of known scalar potentials. 

Considering potential (T4"2l and using its dual shape invariance it is possible to 
discover that for half integer fi V K can be transformed to a direct sum of scalar 
Coulomb potentials. Indeed, its superpartner obtained with using superpotential 
(15"T|) with opposite sign, i. e., = —W^ looks as: 

V& = W^ K + W'^ K + c,= - 1) + k 2 - K (2Ai - l)a 3 ) 1 - (67) 

Considering = ^k,a*+i as ^he main potential and calculating its superpartner 
with using superpotential W^ + i iK we come to equation fIBTj) with fi — > fi — 2, etc. It 
is easy to see that continuing this procedure we obtain on some step the following 
result: 



1(1 + 1) u , 1 



where fl = fi + n, n = —fi — |. Diagonalizing matrix Ci — > a 3 we reduce (I6"5|) to a 
direct sum of attractive and repulsive Coulomb potentials written in radial variables. 
It means that for negative and half integer /i our potential fj4"2]) is isospectral with 
the Coulomb one. 
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In analogous way we can show that potentials ( jH[) with half integer k or integer 
/j, is isospectral with the potential 

V K = X 2 (r (r — 1) sec 2 Xx + 2wtan Asai) , r = - ± /i or r = k, (69) 

which is equivalent to the direct sum of two trigonometric Rosen-Morse potentials. 
Under the same conditions for parameters /z and k potential ( l4"6j) is isospectral with 
the following potential: 

V K = X 2 (r (r - 1) csch 2 (Ax) + 2u coth Xxa 1 ) (70) 

which is equivalent to the direct sum of two Eckart potentials. Finally, potential (146]) 
is isospectral with 



V K = X 2 (r (r — 1) sech 2 Xx + 2wtanh Axcr 3 ) , r = - ± \ hi 2 + - (71) 

provided k is negative half integer. Potential ( 171]) is equivalent to the direct sum of 
two hyperbolic Rosen-Morse potentials. 

Thus for some special values of parameters /i and k we can establish the isospec- 
trality relations of matrix potentials fr4"2"|) -( r4"6"|) with well known scalar potentials. This 
observation is supported by the direct comparison of spectra (I6"2"|) - fl6"2"l) with the spec- 
tra of Schrodinger equation with Coulomb, Rosen-Morse and Eckart potentials which 
can be found, e. g., in [1]. 

Let us note that setting in (H3|) - (|46|) u = we also come to the direct sums of 
shape invariant potentials, namely, Morse, Scraft and generalized Poshl-Teller ones. 
However for nonzero u = and /i, k which do not satisfy conditions imposed to obtain 
f l68|) -f r7T|) the found potentials cannot be transformed to the mentioned directs sums 
using the consequent Darboux transformations. 



9 Ground state solutions 

Let us find the ground state solutions for equations f l57|) with shape invariant po- 
tentials fr42|) - fr45|) . To do this it is necessary to solve equations (158]) where W K ^ 
are superpotentials given in f l34"|) - (l3Tj) . and analogous equation with superpotentials 
f l5^]) -f lo^|) . The corresponding solutions are square integrable two component func- 
tions which we denote as: 



i/j (k,ij,,x) = 




In this section we find the ground state solutions considering consequently all the 
mentioned potentials. 
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9.1 Ground states for systems with potentials (1421) and (1431) 



Let us start with the superpotential defined by equation (1341) . Substituting (134]) and 
(172j) into ( l5"Hj) we obtain the following system: 

^ + (/i - K) ^ + 2^TT^ = ' (73) 
|-^ + " + 1 >! + ^Ti* = - < 74 > 

Solving (1731) for yj, substituting the solution into ([73]) and making the change 
we obtain the equation 

^ +y i- (y2+ ^ )i=0 - (76) 

Its solution is a linear combination of modified Bessel functions: 

i = C 1 K^y) + C 2 I,(y)- (77) 

To obtain a square integrable solution we have to set in (1771) C2 = since 
turns to infinity with x — > 00. Then substituting (1771) into (1751) and using (1741) we 
obtain solutions for system (175j) . (1741) in the following form: 

v , = 2 / k+ % +1 ( 2/ ), e = y K+1 %i(y) (78) 

where y is the variable defined in (T75|) . ux/(2k + 1) > 0. 

Functions (j75]l are square integrable provided parameter k is positive and satisfies 
the following relation: 

k-/z>0. (79) 

If this condition is violated, i. e., 

k - n < (80) 

we cannot find ground state vector using equation (158|) with superpotential (IMI) since 
its solutions ([75]) are not square integrable. But since potential (jl2l) admits the 
dual shape invariance, it is possible to make an alternative factorization of equation 
(j57p using superpotential (j520 and search for normalizable solutions of the following 
equation: 

«> x)i) Q (p,, k, x) = 0. (si) 
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where (and in the following) 

~ a ^ = di +w ^ *t« = -fa+ w ^- ( 82 ) 

Indeed, solving (15T]) we obtain a perfect ground state vector: 

Mv,K,x)=( ( g\ <P = y" + *K\v\(v), £=y^K w _ 1] (y) (83) 

where y = anc ^ v = K ~^ V^- The normalizability conditions for solution (153]) 
are: 

K-fi<l, if k > (84) 

and 

K + > 1, if /t < 0. (85) 
It is important to note that conditions (1791) and (1541) are compatible provided 

k>0, 0<k-//<1. (86) 

Conditions ( 179]) and (185]) are incompatible. 

Thus if parameters /i and k satisfy ( 180]) and ( 154"]) , equation ( I5"7j) admits ground 
state solutions (155]) . If ( 179]) is satisfied but ( 154]) is not true, the ground state solutions 
are given by relations ( 175]) . If condition ( 156"]) is satisfied both solutions ( 175]) and (155]) 
are available. In the special case n = /i + 1/2 solutions ( 175]) and (153]) coincide. 

Notice that our convention that parameter /z is positive excludes the case /z = —1/2 
when potential (T42]) is reduced to a direct sum of Coulomb potentials, see section [8j 

Analogously, considering equation (15^]) with superpotential (155]) and representing 
its solution in the form (172]) with 

£ = 2/^ K £(y), V 5 = y^~ K 0(y) i y = fiexp(-\x) 
we find the following solutions: 

= 2/^m(?/), e = -j/*-"#|„-n(y) (87) 
where v = u/k, + 1/2 and parameters u; and k should satisfy the conditions 

k < 0, /t 2 > u. (88) 

Since potential (143]) does not admit the dual shape invariance, there are no other 
ground state solutions. 
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9.2 Ground states for systems with potentials (1441)— (1461) 



In analogous manner we find solutions of equations fl58|) and flHTj) for the remaining 
superpotentials (13"5"1) - (13"T|) . Let us present them without calculational details. 

Solving equation ([SHD for superpotential (13"6T) we obtain two normalizable solutions, 
the first of which is: 

<Pi = 2/^(1 -y) t ^ B 2Fi{a,b,c; y), 

2u , „ , v (89) 

& = "7^ -n » 1-n » 2 Fi a + l,Hl,c+li!/ . 

K{2jJ, — 1) 

Here 2-^1(0, b, c; y) is the hypergeometric function, 

a = —1—, o = 1—, c = /i ^ U, 

k k 2 

1 7T 7T 



(90) 



and parameters /i and n are constrained by the conditions (179"]) and 

/t + /i>0. (91) 
The second solution is 

<^ 2 = y 1 ^ (1 - y)^2Fx (a, 6, c; y) , 

(2//+1)k /l-y\* 
6 = a^ - V2 (92) 

where variable ?/ is the same as in fl89|) . 

1 - w 7 1 .w 3 , rtfl . 

a = At+7j- 1 - ) b = fi + - + i-, c = //+-, (93) 



and parameters k, /i again should satisfy conditions f J79|) and f l9Tj) . 

Using the dual shape invariance of potential OH]) we can find additional (or alter- 
native) ground state solutions using equation ( 18"T|) with superpotential flo^]) . In this 
way we obtain 

<P\ = Z/ iL ^ ±i (l - y)^2^i (a, 6, c; y) , 

4a; , , i+k+m , , (94) 

5 1 = w^m^T/ ' (1 " » )_,_,Fl (a + L 6 + 11 c + 1: • 

Here variable y and its domain are the same as given in (|90|) . 

2iw , , 

-a = 6 = — , c = l-/c, (95) 
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and parameters fi and k are constrained by the conditions (19~T]) and 

k-h<1. (96) 

The second solution is 

^2 = ^2, 6 = 6 (97) 

where y>2 and £2 are functions defined by equation fl92l) where arguments a, b and c 
of the hypergeometric function differs from fl93|) and have the following form: 



2iw 2iw 1 

a = k , b = k c = k H — , 

2/i + l' 2/i + l 2' 

and parameters k, /i should satisfy conditions ( |9T|) and ( 196]) . 

Thus for potential (|4"4]) we have three versions of constraints for parameters \i and 

k : 

«-//>1, (98) 
k - // < (99) 

or 

0<k-//<1. (100) 

In addition, condition ( 19~T|) should be imposed. 

For the cases fl98l) and fl99l) the ground state solutions are given by equations (189]) , 
([92]) and ([94]), dSZD correspondingly while in the case f lT00|) all solutions flHDl. (192]) . 
([94]) and ([97]) are available. 

For superpotential (157]) we obtain the following solution of equation (15^]) : 

^^(l-^^rV^a,^,!/ 2 ), 

6 = -m + — ^ +1 (i - y 2 ) 1 - K+ ^F 1 (a + 1, b + 1, c + 1; y 2 ) ( ' 01 ; 

c 

where 

a; 1 Xx 

a——, b = a + c, c = - — /i, ?/ = tanh — (102) 

/x A A 

and parameters k, /i, u; satisfy conditions (179]) and f lHHj) . These conditions are com- 
patible iff /i < 0. 

One more solution for superpotential (1371) has components given below: 
& = (1 - y 2 )"~ V^+V^i M, c; y 2 ) , 

f «(2/i- 1)^-1) \ 

<P*={ ^ V^ 2 (103) 

+ — (1 - y 2 ) 1 ""* V"" +2 2^i(« + 1, b + 1, c + 1; y 2 ) 

UJC 

18 



where y is the variable given in ( I102p . 

UJ 3 

a = 1 H , b = c — a, c = /i. (104) 

k 2 

Solution Cl 1 3 p is normalizable provided conditions f l79|) and flHHj) are satisfied. 

Potential ( 1431) possesses the dual shape invariance thus we also should find So- 
lutions of equation (fHTj) with superpotential f[5"4"l) . The explicit expression of these 
solutions can be obtained from (11011) and fllU3j) using the change (j5Dl) . To obtain 
consistent solutions the additional conditions 

/i<0, (2/i + l) 2 > Acu > (105) 

should be imposed instead of (1881) . 

For superpotential (1381) we also have two ground state vectors which solve equation 
(|58|) . The first of them has the following components: 

<Pi = y~^ + ^{l - yY^^ 2 Fi (a,b,c; y) , 

ZllK i (106) 
6 = -- — —y 5 2 + 2^(1 -y)2 2 a» 2 .Fi(a + l,& + l,c + l;y), 

ZUJ + K 

where the parameters and variables are 

a = — i/i, b = i/i, c=-H — 7^0, i/ = -(tanh \x + 1) (107) 

2 k 2 

with k and /i satisfying f l8"Sj) . The second solution looks as: 



ip 2 = y2 2 2«(l-y) 2 a " 2 Fi(a, 6, c; y), 

lu; — / 
2«;/i 



2w - K / 1 - „ , 



2 1 /I,, 2 ,_2 

2/ 1 -f-^(l _j,)f-f-3s 2 F 1 (a + l,6+l,c+l;j/) 



2/i/«(3K — 2u>) 

where y is the variable defined in (I107p . 

lu; lu; 3 u; 

a=----vi, &=--- + ^ c =2"^' (1 ° 9) 

and parameters k, u should satisfy condition (IHHj) . 

Due to the absence of dual shape invariance of potential (jlBl) there are no alter- 
native ground state solutions. 

Thus we find ground state solutions for equation f )57|) with all potentials (j4"2"]) -(l4"5" ]l . 
These solutions are square integrable and correspond to the eigenvalues E K = — c K or 

= — c M where c K and c M are given by equations (j4"Tj) - (j4"9l and fl55]h 
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10 Exited states 



We already know ground state solutions for Schrodinger equations with potentials 
given by relations ( I42p - (|45p . Solutions which correspond to n th energy level can be 
obtained starting, e. g., with the ground state solutions (178|) . (1H7|) . (|89|) . (1921) . fllOip . 
( I103p . f)106p . (I108P and applying equation (l6Tj) . However it is necessary to make sure 
that such defined ground and exited states are square integrable. 

Let us first consider potential ( 142]) which admits the dual shape invariance. This 
invariance enables to make factorization of the corresponding Schrodinger equation 
and find ground state solutions using either superpotential (154"P or (IB^|) . or even both 
of them, depending on given initial values of parameters k and /z. Namely if f J79|) is 
satisfied but (184"]) is not true , the ground state solutions are given by equations flT2"]) 
and f JTHj) and exited states are given by equation f l61~]) . 

The corresponding energy levels are given by equations f[6"2"]) and (ISI)]) . 

If parameters /i and n satisfy (ISU]) and one of conditions (1541) or (1541) . equation 
( |57p admits ground state solutions (183]) and exited states are defined by the following 
equation: 

$ n {K,ii,x) =a+ M a+ +1)M ---a+ +n _ 1)M '0o(K + n,/i,a;) (110) 

where a+ = — ^ + and W^ K is the alternative superpotential given by fl52|) . 
The corresponding energy levels are given by formulae (162]) and (166"P 

If condition f )86|) is satisfied both versions of solutions and energy levels given 
above are available. In the special case n = fi + 1/2 solutions f!78|) and fl83l) coincide. 

Let us start with the ground state solution f[75|) . Its normalizability is almost 
evident since the modified Bessel function has the only singular point, namely, y = 0, 
and decreases exponentially at infinity. Moreover, at y — there is the inverse power 
singularity, i. e., K u ~ \ with u = jioiu = ji + l which is perfectly compensated by 
the multiplier y K+1 provided k satisfies fTTUl) . 

Calculating the first exited state (IBU]) we can use relation flo^]) where « — y k + 1, 
thus 

Mk,»,y) = (W K+1 ^ + W K;tl )ip (K + i,v,y). (m) 

Again we recognize a good behavior at the singularity point y = since (W^+i^ + 

W KjfJ ) ~ J H and ^o(« + 1, A*, y) ~ #o(«, /i, 2/)- 

Let us suppose that the wave function corresponding to n th exited state is regular 
at y — and has the following form 

^n(K,/J.,y) = Wt/j (K + n,fi,y) (112) 

where ipo(K + n, /i, y) is the ground state solution given by relations f T72l) . fTTBl) and 
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is a matrix depending on y and k. Then 



'(p n+ i(K,n,y) = (—^- + W K ) 4> n (K + n + l,fi,y) 

V " J x (H3) 
dW \ 

W K W + WW K+n+l + n + 1, /i, y) 



dy 



where we use the fact that i/)q(k + n + 1, /i, y) solves the equation (158]) where k — >• 
/t + n + 1. 

Using ( |113p it is not difficult to show that if function tp n (n + l,n,y) be square 
integrable then ip n+ i(n, /z, y) is square integrable too. Indeed, at the neighborhood 
of the singularity point y = the ground state functions are related as ipo(K + n + 

l,/x,y)~yVo(« + n,/i,!/),and + W K W + WW K+1 ^ ~ \W . Thus ^ n+1 ( K , /i, y) 

is regular at y = provided ^(^j Z 1 ; 1/) be regular. Since for n = 0, 1 our supposition is 
fulfilled we conclude by induction that wave functions i/j n (n, /i, y) (l61i) are normalizable 
for any n. 

By direct repeating the above speculations we can prove the square integrability of 
ground state vector (1531 and exited state vectors given by relation (11 lop . In fact the 
only thing we need is to change ip n (n, [i, y) and W Kjll by their counterparts ip n (n, /i, y) 
and W w 

In complete analogy with the above one can prove the normalizability of the exited 
states for the case when the superpotential is given by equation (j3"5|) . However there 
is an essentially new point which is generated by condition (|88p . The think is that 
solutions (18"?]) being well defined for k and /i satisfying conditions (1791) and (1881) . can 
loose their square integrability after the change k — > k + n for a sufficiently large n. 
Namely, in order to obtain a normalizable solution ipo{ K + n , %) we have to ask for 
(k, + n) 2 > u > 0. Since k is negative we have the following restriction for n: 

ti<\k\- ^fbj. (114) 

It is possible to show that if ipo(n + n,fi,x) is not normalizable the same is true for 
exited states (15Tp . 

Let us consider the ground state solutions (18"9"|) . This solution like all the re- 
maining solutions ([92]), (JM]), (J9JJ), ffTOT]) . ffT03]) . dTM]) . ( 1108]) is expressed via linear 
combinations of the following elements: 

y A (l-y) B 2 F 1 (a,b,c;y) (115) 

where parameters a, b and c are given by equation (190]) . A = ^p,.B = for 
component etc. 

In accordance with its definition, variable y belongs to the interval [0, 1] and there 
are two points which are "suspicious w. r. t. singularity", namely, y = and y — 1. 
In order the solution to be regular (and equal to zero) in these points it is necessary 
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and sufficient to ask for A > 0, B > and dte(B + c — a — b) > 0. Exactly these 
conditions generate restrictions (|88|) for parameters k which guarantee the solution 
normalizability. The same is true for solutions (J92]), fTT06]) . f lT08|) . ffToB and ffT03]) . 

To analyze solutions for exited states we rewrite superpotential (136|) in terms of 
variable y: 

W K ^ = X («(22/ - 1) + ^(1-?/)^ + ^<x 3 ) . (116) 

We see that W K is nonsingular at y = and y = 1, the same is true for W K+n ^ for 
any natural number n. Functions (l89il are still regular at these points if we change 
k K + n, thus we can again apply relations (Illip - (lll3p to prove the normalizability 
of the corresponding solutions (161 p for arbitrary n. 

In a similar way we can prove the square integrability for solutions (1611) cor- 
responding to ground state solutions ([92}, flU}, ([97]), ffTOTD . ffT03l) . ffT06]) and ffT08|) . 
However in these cases we again have constraint ([55]) which generates restriction (11141) 
for the number n enumerating the exited states f[6T]) . Analogously, starting with con- 
dition (I105p we come to conclusion that solutions (IllOp for the Schrodinger equation 
with potential (j4"o"|) are square integrable provided quantum number n satisfies the 
condition 

n<H-v^-^ (117) 

Thus for any fbced k, /i and u equation f [57l) with potentials f[43]) f[45]) and f[46]) 
describes a system which has a finite number of states with discrete spectrum. These 
states are enumerated by non-negative natural numbers n satisfying condition (I114p 
ordHH). 

The systems with potentials fj4^]) - fj4T]) can also have states with continuous spec- 
trum. In particular, such states should change the bound states when conditions 
(11141) and (I117p are violated. Analysis of the states with continuous spectra lies out 
of frames of the present paper. 

11 Discussion 

Generalizing the supersymmetric PS problem we find a family of matrix potentials 
for Shrodinger equation satisfying the shape invariance condition. In this way we 
find five exactly solvable problems for systems of coupled Shrodinger equations. The 
related matrix potentials are given by equations f ]42]) - (]46]). 

Let us stress that we present the completed classification of shape invariant super- 
potentials of the generic form ([9]) where P and R are hermitian matrices of arbitrary 
finite dimension and Q is proportional to the unit matrix. Namely, we show that such 
objects can be reduced to direct sums of known scalar superpotentials and superpo- 
tentials presented in section [5j 
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The found potentials include parameters A, k, /z and u whose possible values are 
restricted but quite arbitrary. Moreover, parameters u in fH2|) and /i in fj43l) can be 
reduced to unity by scaling and shifting the independent variable x correspondingly. 

Taking into account all possibilities enumerated in f!79|) . (184|) . fl85|) and fl98|) - fll00l) 
we conclude that in the case of potentials PSl) . (jUJ) and (S5J) there are discrete 
spectrum states for all real values of arbitrary parameters A, k and /i except the case 
k = /i. Parameter u; can be constrained by equations (188]) or (llOSp . 

Potential (H2l) is a slightly generalized effective potential for the PS problem. 
Moreover, these potentials coincide for a particular value /z = of arbitrary parameter 
fi. However, if fi ^ potential (I42p is not equivalent to the potential appearing in 
the PS problem and corresponds to a more general interaction in the initial three- 
dimension problem. 

At the best of our knowledge the remaining potentials fH3|) - fH6l) are new. The 
related Schrodinger equations can be integrated using tools of the SUSY quantum 
mechanics. The corresponding spectrum and eigenvectors are given by equations 
flo^]) - fl^oT) and (RJTj) or fiil 01) while the ground state solutions are discussed in section 
M solutions Notice that the "matrbc supersymmetry" has a new feature in comparison 
with the Standard (i. e., scalar) one. Namely, matrix models with shape invariance can 
have degenerated ground states in spite of that there exists a normalizable solution 
for equation (1581) . Example of models with such specific spontaneously broken SUSY 
is given by the Schrodinger equation (loTI) with potentials (S3D (SSJ) • 

Mathematically, there are natural reasons for appearance of a zero energy doublet 
of the ground states in systems with the matrbc supersymmetry. The thing is that 
equation fl58|) is a system of two the first order equations whose solutions are linear 
combinations of two functions while in the ordinary SUSY quantum mechanics we 
have a one first order equation for ground states. For potentials fT4"2"l) and (T4"3"l) only 
one of these functions is normalizable but for potentials (r4"4"l) - (T4"3"l) there are two ground 
state solutions. 

Let us note that existence of zero energy doublets of the ground states was already 
registered in periodic quantum systems, see [18] and [19] for discussion of this phe- 
nomenon. In this connection it seems to be interesting to extend our approach to the 
case of periodic systems. Formally speaking, the only new constructive elements of 
potentials fH2l) - fj46l) in comparison with the Standard scalar shape invariant potentials 
are matrices a± and 03 which are involutions anticommuting between themselves. In 
fact the nature of these involutions is not essential for deducing the shape invariant 
potentials, and many of the results discussed in present paper can be generalized to 
the case of another involutions. For example, it is possible to change the mentioned 
matrices by reflection and shift operators which also can be anticommuting involu- 
tions being applied to functions with an appropriate parity and periodicity. In this 
way it seems to be possible to extend the list of potentials which admit supersym- 
metry including shifts of arguments [20]. A classification of anticommuting discrete 
symmetries and the corresponding supersymmetric versions of the Schrodinger and 
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Pauli equations can be found in [22] . 

An interesting phenomena which appears to be typical for systems with matrix 
SUSY is the dual shape invariance discussed in Section [6j It enables to impose 
much less restrictive constrains on parameters of potentials then the ordinary shape 
invariance. In addition, it can be used to explain the insensibility of the spectra 
( I62p -( 1fi2"|) . (I65p on parameter /i. Namely, hamiltonians with shifted fi should be 
almost isospectral thanks to the dual shape invariance, which is incompatible with /i- 
dependence of energy values (excluding the exotic case when these values are periodic 
functions of /i). 

For some values of parameters [i and k the additional branch of spectrum caused 
by the dual shape invariance can appear. In particular it is true for potential (|42j) 
with /i = and < k < 1/2. Enhanced analysis of such potentials was made in 
paper [12] . In the present paper we slightly refine results of [12J. 

Let us note that the dual shape invariance can be recognized for two potentials of 
the ordinary SUSY quantum mechanics, namely, for the trigonometric Scraft 1 and 
generalized Poshl-Teller potentials: 

V\ = (k(k — 1) + /i 2 ) sec 2 x + — 2k) sec Xx tan Xx, 

V2 = (k(k — 1) + fj 2 ) csch 2 (x) + — 2k) coth Xx csch x 

both of which admit symmetries (150]) . The corresponding energy spectra is /z-independent 
also, and the dual shape invariance can be used to explain this phenomena. 

It is shown in section [S] that for some values of parameters /1 and k the matrix 
potentials (I42 p — (j46p are isospectral with direct sums of one dimensional shape invari- 
ant potentials. Unfortunately, in this way we cannot establish isospectrality with the 
reflectionless hyperbolic Poschl-Teller (HPT) system which has a lot of interesting ap- 
plications and admits a hidden (bosonized) nonlinear supersymmetry [21] . A natural 
question arises if there are other matrbc superpotentials which can add the list given 
by equations (J54"]) - (}H7|) and may include a matrix counterpart of the reflectionless 
HPT potential? 

In sections 3-5 we show that irreducible matrbc potentials of generic form fl9]) are 
exhausted by known scalar ones and 2x2 matrbc operators given by equations ( )34|) - 
( 13 7p . Of course we restrict ourselves to the shape invariance of type 1 when variable 
parameter k is changing by shifts. 

Nevertheless it is possible to search for matrbc superpotentials in a more general 
approach, when matrbc Q in fl9]) is not restricted to be proportional to the unit one. 
Let us present an example of such superpotential: 



where c is a a constant. The corresponding potential V K ()6]) is shape invariant and 
has the following form: 




V K = {Ak 2 - 1) 



x 2 + c 2 — 2cxa 3 



+ 



4(x 2 - c 2 ) 2 



c 2 — x' 
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while the energy spectrum is given by equation (16"2"1) . The other examples (which, 
however, are restricted to a linear dependence of W K)fi on variable parameter k) can 
be found in paper [TTj . 

The problem of classification of matrix potentials fl9]) with generic hermitian ma- 
trices Q, P and i? is a subject of our contemporary research. 
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